Influence of nonlocal electrodynamics on the anisotropic vortex pinning in YN12B2C 
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We have studied the pinning force density F p of Y N12B2C superconductors for various field orien- 
tations. We observe anisotropies both between the c-axis and the basal plane and within the plane, 
that cannot be explained by usual mass anisotropy. For magnetic field H || c, the reorientation 
structural transition in the vortex lattice due to nonlocality, which occurs at a field H\ ~ IkOe, 
manifests itself as a kink in F P (H). When H_Lc, F p is much larger and has a quite different H 
dependence, indicating that other pinning mechanisms are present. In this case the signature of 
nonlocal effects is the presence of a fourfold periodicity of F p within the basal plane. 



After Kogan and co-workersEI showed in 1996 that non- 
local electrodynamics effects are relevant in high-K su- 
perconductors in a broad H-T domain, several works 
demonstrated that nonlocality give rise to unusual prop- 
erties of the vortex matter. Those studies were done 
mainly on the family of compounds RN12B2C ', where 
R = Lu, Y, Tm, Er, Ho, Dy. These borocarbides are 
very well suited to study nonlocal effects due to their in- 
termediate k values (~ 10 - 20) that result in a wide field 
range where the extended London description holds, their 
high transition temperature T c and the possibility to fab- 
ricate clean single crystals with large electronic mean free 
path. Deviations of the equilibrium magnetization from 
the local London prediction] 2 ] structural phase transitions 
separating a variety of exotic rhombic and square vortex 
lattices (FLL)£TLL3 and fourfold anisptcapies in the basal 
plane of these tetragonal materialsUta provide strong 
support for the nonlocal scenario. 

Although the consequences of nonlocality on the 
equilibrium properties of the superconducting vortex 
matter is by now convincingly established, very little is 
known about its effects on the nonequilibrium vortex re- 
sponse^ For both YNi 2 B 2 C and LuNi 2 B 2 C Eskildsen 
et al.E 3 ] have shown that, for H || c-axis and well above 
the field H 2 ~ IkOe marking a transition from a rhombic 
to a square vortex lattice, the pinning properties sharply 
disagree with the expectations of the Larkin-Ovchinnikov 
collective pinning theory for a triangular lattice. The au- 
thors associated the discrepancy to the different (and to 
large extent unknown) elastic properties of the square lat- 
tice. However, neither the irreversible response in other 
field orientations, where the lattice geometry is different, 
nor the influence of the lattice transitions themselves on 
vortex pinning have been experimentally explored until 
now. It is clear that vortex pinning, which involves dis- 
tortions from equilibrium vortex configurations, should 
be affected by the symmetry changes in the lattice. 

In this work we report a study of nonlocal effects on the 



vortex pinning of a single crystal of Y Ni 2 B 2 C . For H || 
c-axis the field dependence of the pinning force density 
F p exhibits a distinct kink at a field H* which coincides 
with the first order reorientation transition in the vortex 
lattice. The angular dependence of H* (for H inclined up 
to 30° from c) provide further support for the association 
of the kink with the FLL transition. For H in the ab- 
plane, F p is an order of magnitude larger and the field 
dependence is quite different. We discuss the possible 
origin of this anisotropy. Finally, we observe fourfold 
oscillations in the angular dependence of F p in the basal 
plane which also arise from nonlocal effects. 

The single crystal measured in the present study, with 
a mass of 17 mg, dimensions ~ 2 x 2.5 x 0.5mm, and 
critical temperature T C =14.5K, is the same one used in 
previous worksElid. Electrical resistivity data yielded a 
residual resistance ratio of 10 and an estimate of the elec- 
tronic mean free path I ~ 300A Magnetic measurements 
were carried out in two Quantum Design SQUID magne- 
tometers. All data were taken in a MPMS-7 (H < IT) 
except the angular dependence in the basal plane, which 
was studied in a MPMS-5S (H < 5T) using a home made 
rotating system. til In the superconducting mixed state, 
isothermal magnetization loops were recorded. The hys- 
teresis width AM (if) was used to calculate the critical 
current density J c from the Bean's critical state model. 

Figure 1(a) shows the pinning force density F p = J c x 
B as a function of H || c-axis at several temperatures. 
At T = 3K, after an initial increase at very low fields 
(not shown), F p (H) reaches a maximum at a remarkably 
low reduced field (H/H c2 < 0.01), and then decreases 
approximately linearly in log(H). At a field H* ~ IkOe 
a clear kink is observed. Above this field F p is still linear 
in log(H), but the slope dF p / dlogH has a smaller value. 
At a higher field Hf, ~ lOkOe a small bump is visible. As 
discuss below, we associate Hb with the peak observed by 
Eskildsen et al£j in LuNi 2 B 2 C. Finally, close to H c2 we 
observe the usual peak effect attributed to the softening 
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of the flux line lattice. n 

In the framework of Kogan's modelu two structural 
transitions in the vortex lattice occur in YNi 2 B 2 C, 
namely a first order reorientation transition between two 
rhombic lattices at a field ti\ and a second order tran- 
sition from the rhombic to a square lattice at a higher 
field |££p. According to a recent analysis by Knigavko 
et al.Ej the reorientation at Hi really consists of two 
second order transitions, but they take place in such a 
narrow field range that may be hard to resolve them 
experimentally. Small angle neutron scattering (SANS) 
experiments^ (performed on a crystal of the same batch 
as ours) confirmed Kogan's predictions. For H || c, a 
jump in the apical angle (3 of the rhombic vortex lattice 
(discontinuous within the resolution) occurs at Hi ~ 1 
to 1.25/cOe, and the lattice becomes square {j3 = 90°) at 
H 2 ~ 1.25 to 1.5kOe. The coincidence of our H* with 
this field range make it tempting to associate the jump in 
dF p /dlogH with one of these transitions. However, the 
small difference between Hi and H2 make it difficult to 
determine which one of them generates the kink. 

To solve this question we need to modify the experi- 
mental conditions in such a way that Hi and H2 change 
and split apart. The simplest way is to change T. Fig. 
1(a) shows that as T increases H* remains almost con- 
stant, while the slope change progressively washes out. 
This last observation is consistent with the association 
of H* with nonlocality, as those effects tend to disap- 
pear as T approaches T c . On the other hand, theoretical 
expect atipnso as well as some qualitative experimental 
evidenced suggests that H 2 increases with T well below 
T c , in contrast with the observed behavior of H*. This 
seems to indicate that the kink is not related with H 2 . 

To our knowledge there are no experimental data on 
the temperature dependence of Hi . The only theoretica.1 
hint comes from the recent study of Knigavko et al.E-3. 
It is shown there that Hi in units of the scaling field 
<I>o/(27rA) 2 is a decreasing function of the strength of the 
nonlocal effects (which is parametrized by two averages 
over the Fermi surface, n and d), and consequently must 
increase with T. Disregarding a weak temperature de- 
pendence in k, the scaling field is proportional to H c2 , 
which is easily accessed experimentally. This suggests to 
analyze F p as a function of the reduced field h = H/H c2 . 

In Figure 2 we replotted the data of Fig. 1(a) as 
f p = F p /FP eak vs h. Here F^ eak is the maximum of 
Fp at the peak close to H C 2, which occurs at the same 
reduced field h p — 0.8 for all T. The overlap of the data 
at different T in the h range of the peak-effect means 
that the pinning mechanism responsible for it satisfies a 
scaling law F p {H,T) = FP eak (T)f p (h), as seen in many 
systemsJlj On the other hand, the fact that at lower h 
the fp(h) data at various T do not collapse on a single 
curve indicates that the dominant pinning mechanism(s) 
in that range is(are) different from that originating the 
peak-effect. As already mentioned, we identify the bump 
at intermediate field with the peak reported by Eskild- 



sen et al.Ej The fact that this bump occurs at the same 
hj, ~ 0.18 for all T is consistent with the interpretation of 
this feature as resulting from the stiffening of the vortex 
matter due to the increase of the shear modulus Cee^AS 
this elastic constant depends on field only through hE3 

In contrast to h p and hi,, Fig. 2 shows that h* 
II" II,'} is temperature dependent. This lack of scaling is 
expected for any feature arising from nonlocal effects, as 
they involve a new characteristic field Hq oc n~ 2 , whose 
T dependence is quite different from that of H c2 . In par- 
ticular, the fact that h* increases with T is consistent 
with the theoretical prediction for hi = Hi/H c2 as dis- 
cussed above. So, from the temperature dependence of 
H* we conclude that it is probably associated with Hi. 

Further evidence that H* is a signature of the reori- 
entation transition at Hi arises from the angular depen- 
dence. SANS results showo that H 2 is very sensitive to 
the field orientation, while Hi is much less so. When H 
is tilted by an angle O = 10° from the c-axis, Hi re- 
mains unchanged within the resolution, while H 2 shifts 
up significantly. Extrapolation of the data (the maxi- 
mum measured field was 2.5/cOe) gives H 2 ~ 2.7kOe. 
For O = 30° only four SANS data points are available. 
Now (3 increases only weakly with H , and extrapolation 
suggests that H 2 is well above AkOe. The inset of Fig. 2 
shows F p vs log(H) at T = 3 and 5K, for O = 0°, 10° 
and 30°. The 10° curves are almost identical to H || c 
and the kink is clearly visible at the same H*, indicating 
that it is related to Hi rather than to H 2 . In the O = 30° 
curves H* is slightly shifted up, while according to the 
SANS data Hi has either disappeared or shifted down to 
H < IkOe. However the SANS evidence in this case is 
inconclusive, as it is based on a single data point and it 
may be affected by phase coexistence. 

The very low J c for H || c indicates that pinning corr,e-| 
lation volumes are large, as confirmed by SANS results.^ 
Thus, the elastic properties of the vortex lattice must 
play a key role in the pinning. While the tilt modulus 
C44 = BH/Att should be rather independent of the lat- 
tice details, it is easy to see that Cqq depends on (3. As 
(3 overgoes a discontinuous jump at Hi, it is natural to 
ascribe-the kink to an abrupt change in 6*66 ■ It has been 
arguedEj that below Hb the collective pinning description 
is valid in YNi 2 B 2 C. In the simplest scenario this im- 
plies that F p oc C^Cqq. Using this expression and the 
experimental F p {H) we computed C§§{H), and indeed 
found a kink at H* as expected. The surprising result 
is that Cqq decreases monotonically with H both below 
and above H* , in contradiction with the standard behav- 
ior. It must be noted, however, that the shear properties 
for (3 7^ 60° are anisotropic, thus Cqq has more than one 
component and the expression for F p should be modified. 

We now turn to the pinning properties for HJ_c. Figure 
1(b) shows Fp(H) for H || [100] and H || [110] at various 
T. The behavior in these orientations is very different 
from the H || c case. At T = 3K a broad maximum 
occurs at fields H max ~ 8kOe for H || [100] and slightly 
less for H || [110]. Over most of the field range F p is 
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much larger than for H || c, the maximum anisotropy 
being a factor of ~ 20 at H rnax . As T increases, H max 
decreases while the out-of-plane anisotropy in F p remains 
approximately constant. Another feature visible in Fig. 
1(b) is the in-plane anisotropy, F p [100] > F p [110]. The 
difference is small, but clear and systematic, and holds 
at all H below the peak effect. This in-plane anisotropy 
smoothly decreases with T. 

The large difference in F p between the c-axis and 
the basal plane is surprising, considering that the mass 
anisotropyc3 is very small (< 10%) as confirmed by the 
H C 2 values in Figs. 1(a) and (b). This rules out ex- 
planations based on the anisotropic scaling frequently 
used in high T c superconductorsHI And of course mass 
anisotropy cannot account for any angular dependence 
in the basal plane of this tetragonal compound. We will 
now discuss the possible sources of both the out-of-plane 
and in-plane anisotropy. 

From the analysis of the normal state magnetization 
M ns we foundtil that this crystal contains a small amount 
of magnetic impurities, about ~ 0.1% of a magnetic rare 
earth substituting for Y. The moments lay in the basal 
plane and produce a Curie tail that, for values of the 
scaling variable H/T up to ~ IkOe/K, is isotropic within 
the plane. For higher H/T a fourfold in-plane anisotropy 
develops. This is shown in the inset of Fig. 3, where 
5M ns = M ns [110] — M ns [100] is plotted as a function of 
H/T. These localized moments are a potential source of 
anisotropic vortex pinning. However, although we cannot 
completely rule out this possibility, several observations 
suggest that this is not the case. 

In the first place, as the alignment of the moments 
increases with H, the directional pinning should become 
more efficient as the field increases and thus the out- 
of-plane anisotropy in F p should increase monotonically 
with H. This is not in agreement with the observations. 
At T = 3K the anisotropy factor is already large (~ 6) 
at H <~ IkOe, then it grows with H, maximizes at H ~ 
H m ax and decreases again for higher fields. In particular, 
F p is almost isotropic (within ~ 10%) in the peak effect. 

Another clue comes from the in-plane anisotropy. As 
M ns is isotropic within the plane for H/T < IkOe/K (see 
inset of Fig. 3), we would expect F p to be also isotropic 
in the plane below that threshold. This is again contrary 
to the results: Fig 1(b) shows that at all temperatures 
3K < T < 9K the in-plane anisotropy in F p is already 
visible at fields well below the H/T threshold. 

The lower panel of Figure 3 shows the in-plane an- 
gular dependence of the pinning force density, F p (ip), at 
T = 7K for several H, and the upper panel shows M ns (ip) 
for H = 45fcOe. This figure clearly demonstrates that F p 
exhibits a fourfold anisotropy for all H > 2kOe, while at 
this temperature the threshold for the anisotropy in M ns 
is ~ 7kOe. The behavior of F p (ip) is rather complex: 
in addition to the main peaks at [100] and [010], sec- 
ondary maxima are visible at [110] and (110). Further 
evidence that the localized moments are not the source 
of anisotropic pinning is that the fourfold oscillations in 



M ns and F p have opposite sign, as is evident by compar- 
ing the lower and upper panels. The presence of four- and 
higher-fold oscillations in Fig. 3, as well as the nontrivial 
T and H dependencies in Fig. 1(b), rule out the possi- 
bility that the in-plane anisotropy could be an artifact 
due to a slight misalignment between H and the basal 
plane. Such misalignment would only produce two-fold 
oscillations, which are not observed. 

Once the possibility that the in-plane anisotropy in F p 
is due to the magnetic impurities has been ruled out, we 
must consider other alternatives. While it is very unlikely 
that any type of crystallographic defects could account 
for the complex in-plane variations in F p , nonlocal effects 
can provide a natural explanation for them. Due to non- 
locality, the geometry of the vortex lattice depends on 
the orientation within the plane. We can again argue (as 
in the H || c case) that Cqq depends on the lattice geom- 
etry, and that such dependence must be reflected in F p . 
Interestingly F p has local maxima at the high symmetry 
crystallographic orientations [100] and [110], where Cqq 
is expected to exhibit local minima. The fact that the in- 
plane anisotropy disappears as T approaches T c provide 
further support to this interpretation. 

The origin of the out-of-plane anisotropy in F p is less 
clear. On one hand, it seems too large to be ascribed to 
nonlocal effects. On the other hand, a simpler explana- 
tion for it could be the presence of some still unidenti- 
fied anisotropic pinning centers, such as planar defects. 
Clearly, the unexpected complexity of the angular de- 
pendence of F p involves other mechanisms besides non- 
locality, and further studies are necessary to sort out the 
various sources of anisotropic pinning in this material. 

In conclusion, F p in YNi 2 B 2 C exhibits a rich 
anisotropic behavior that sharply contrasts with its small 
mass anisotropy. Nonlocal electrodynamics influences 
pinning via the unusual behavior of the shear modulus in 
non-hexagonal lattices, cither through continuous varia- 
tions (fourfold basal plane anisotropy) or abrupt jumps 
(the kink at Hi ) . A complete understanding of the vor- 
tex irreversible response will require a deeper knowledge 
of the elastic properties of non-hexagonal lattices and the 
extension of the pinning models to those structures. The 
origin of the large out-of plane pinning anisotropy, which 
almost certainly does not originate on nonlocality, also 
deserves further investigation. 
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Fig.l: Pinning force density F p vs H for YNi^B^C 
single crystal with (a) H || c-axis and (b) H || [100] 
(open symbols) and H || [110] (solid symbols) at several 
temperatures. 

Fig. 2: Reduced pinning force density F p /F p lax vs 
H/H C 2 at several temperatures with H || c-axis. The 
inset shows F p versus H for H || c (G = 0°) and tilted 
10° and 30° from c. 

Fig. 3: Upper panel: the normal state magnetization 
M ns at T = IK and H — 45kOe, as a function of the 
angle ip between the applied held H (contained in the 
ab basal plane) and the a axis. Lower panel: pinning 
force density F p (ip) at T — 7K and several fields. In- 
set: oscillation amplitude of the in-plane normal state 
magnetization, SM ns = M ns [110] - M„ s [100], vs H/T. 
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